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Studying neural coding through neural energy is a novel view. In this paper, based 
on previously proposed single neuron model, the correlation between the energy 
consumption and the parameters of the cortex networks (amount of neurons, coupling 
strength, and transform delay) under an oscillational condition were researched. We found 
that energy distribution varies orderly as these parameters change, and it is closely 
related to the synchronous oscillation of the neural network. Besides, we compared 
this method with traditional method of relative coefficient, which shows energy method 
works equal to or better than the traditional one. It is novel that the synchronous activity 
and neural network parameters could be researched by assessing energy distribution 
and consumption. Therefore, the conclusion of this paper will refine the framework of 
neural coding theory and contribute to our understanding of the coding mechanism of 
the cerebral cortex. It provides a strong theoretical foundation of a novel neural coding 
theory — energy coding. 
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INTRODUCTION 

Studying neural coding and decoding is one of the most impor- 
tant, dynamic, and challenging fields of neuroscience (Amari 
and Nakahara, 2005; Gazzaniga et al., 2009). Yet, there are 
still large obstacles in researching this problem (Purushothaman 
and Bradley, 2005; Natarajan et al, 2008; McLaughlin, 2009) 
because there is no effective theory of neural coding and decod- 
ing to investigate global brain activity (Laughlin and Sejnowski, 
2003). The practical resources required to solve this problem 
are insufficient. Apart from technical limitations of experimen- 
tal neuroscience, most neural computation is confined to several 
conventional modes based on H-H equations. In the view of neu- 
roscience history, valid theories, and proper assumptions usually 
guide research to elucidate new phenomena (Abbott, 2008). 

The neural activity and operations of the brain involve the 
principle of depleting minimized energy and maximizing sig- 
nal transmission efficiency (Laughlin and Sejnowski, 2003). It 
is logical to exploit these operations to survey the neural cod- 
ing and cognitive functions of the brain. Indeed, a series of 
achievements have been obtained by research that uses energy 
theory to explore the correlation between nervous energy and 
cognitive behavior (Levy and Baxter, 1996; Raichle and Gusnard, 
2002; Wang and Zhang, 2006, 2007, 2012; Wang et al, 2008, 
2009). The most remarkable finding is a new functional mech- 
anism of neurons. That is, while firing action potentials, neurons 
firstly absorb energy, namely oxyhemoglobin, and then con- 
sume energy, namely deoxyhemoglobin. This mechanism chal- 
lenges the traditional energy theory that neurons only consume 
energy. We found that: (1) action potentials strictly corresponds 
to nervous energy, and the strict correspondence between neu- 
ral network activity oscillation and the exiting energy flow. 
This provides a novel method to investigate neural information 



processing through nervous energy, which we term neural energy 
(Wang and Zhang, 2006); (2) The nerve energy encoding can 
automatically assign the nervous energy corresponding to each 
frequency of synchronous oscillations in accordance with the 
allocation ratio in the total energy of the network and sepa- 
rate each neuronal group from the network. This property is 
significant for functional neural network analysis, especially for 
higher cognitive functions (Wang and Zhang, 2012); (3) The 
neural energy can be superimposed, which greatly facilitates 
the analysis of complex neural networks with large dimensions 
and strong non-linearity. Moreover, the superimposed energy 
can reflect the state of synchronous oscillation at the current 
moment. Synchronous oscillation of cortical processing is closely 
related to the formation of conscious cognition (Wang and 
Zhang, 2011); (4) This neural function mechanism can illustrate 
the biological phenomena that have remained indecipherable so 
far (Wang and Zhang, 2012), and theoretically demonstrated 
that neural energy is indeed strongly coupled with variation in 
blood flow (Moore and Cao, 2008; Lin et al, 2010), which, 
consequently, provides a theoretical basis for further neural 
modeling. 

The cognitive and conscious activity of brain actually func- 
tions through control of the transition between structural 
networks and cognitive networks. It is necessary to study a sys- 
tematical method which addresses how the neural computing 
of local network relate to the overall mode of neural activity. 
Synchronous oscillation which caused by dynamic equilibrium 
among different regions and functional combination originates 
from polyrhythmic motions of neural populations in encephalic 
regions. This paper is preliminary and basic prob into cognitive 
neurodynamical mechanism of synchronous oscillation from a 
novel aspect. 
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Yet, this method it is confined to the structural connection of 
neural networks rather than functional connection in research 
with the energy method (Wang and Wang, under review). 
Though it is an important issue, there are no existing reports 
regarding neural energy coding in the cerebral cortex and how 
structural connections convert to functional connection. The aim 
of this paper is to investigate the energy- coding mode of struc- 
tural networks by adjusting the networks' parameters based on 
structural connections. It may provide strict theoretical bases for 
further research about networks' energy derivation and their cor- 
responding coding mode during the conversion from structural 
to functional connections. 

Nervous energy is related with neural signal coding (Wang and 
Zhang, 2006, 2007; Wang and Wang, under review), and net- 
work behavior is also inherently correlated to the parameters. 
Therefore, the total energy caused by the overall network activ- 
ity is interwoven with its parameters (Dhamala et al., 2004). In 
other words, network activity and behavior controlled by param- 
eters (Ghosh et al, 2008), e.g., neuron number, coupling strength, 
and the transform delay, can be described by the network energy 
as it is in oscillation. 

Previously, we showed that as a neuron fires which occur the 
action potential, it firstly discharges the energy in storage (nega- 
tive energy), and then oxyhemoglobin supplies additional energy 
(Wang and Zhang, 2012; Wang and Wang, under review). The 
characteristic about negative energy can effectively reflect neu- 
ronal synchronization. In this paper, this property is taken further 
to investigate the correlation between network parameters and an 
energy- related index, which is the ratio of discharged energy to 
accumulatively consumed energy. 

BIOPHYSICAL MODEL OF NEURONS AND THEIR ENERGY 
FUNCTION 

In order to simulate the energy neuron network coding, we 
propose a novel biophysical model of a neuron comprised of 
voltage and current sources and an inductor (Figure 1). The volt- 
age source is the difference in concentrations of various extra- 
and intracellular ions that induce ion movement. The current 
source is formed by ion concentration gradients and the stim- 
ulated by neighboring neurons. Moreover, when charged ions 
(iC + ,N~+, C+ 2 ) flow into and out of the ion channels, a self- 
induced loop current is formed, which is equivalent to a voltage 





























i i_ 







FIGURE 1 | Biophysical model of a neuron. 



source. The current and the voltage sources can cause energy loss, 
r m and v§ m describe the resistance across P = 5Zm = l ^ m 

and 17, 

respectively. The current and voltage sources are not affected in 
the same position, and film resistor can be divided by r\ m , r2 m , 
and f3 m . The sum of the input current from connected neurons 
describes the coupling relation: 

n 

Im = km + XI [* 0m ^ ^ ^famO' f j- l))] 

7=1 

+iom(tt) sin(co m (n)(f - t n )) (1) 

Where i m \ is the current to maintain the resting membrane 
potential, z'om is the total effect generated by the current stimu- 
lation of peripheral neurons, and oo m is firing frequency. 
The circuit equation corresponds to Figure 1: 



U m = TQ m lQ m + T\ m l\ m + L m I\ m 

lom — I\m ~ Im ~~r~ ~J M ~~r~ C m L7rjm (2) 



and 



Lm-^lrn ~~r~ Tlm^m — K\ m Uo m + K2m ^0m ^2m^m ? (3) 



here K lm = C m (r 2m + r 3m + k 2m = 1 + ^ 

The total power consumed by the neural network composed of N 

neurons is 



p=J2 p » 



(4) 



m = 1 



where the power consumed by the rath neuron is 

P m = d\ m \J^ m + d,2m Uom H~ d^ m UQ m UQ m + dA tm U^ m 
+ d$ m U() m + d(y m 



(5) 



here d\ m , dim, ^3m> <^4m> 4m) and d^ m are the same as in Wang 
and Zhang (2011). 

According to the theory of circuit, energy is the integral of 
power. All the "power" in our subsequent simulations are instan- 
taneous power. In this sense, the power of neuron and the energy 
of neuron is consistent concepts. 

According to Equations (2)-(5), we obtain 
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The simulation result of membrane potential and power con- 
sumed during firing is shown in Figure 2. This novel biophys- 
ical model can simulate neuronal firing well. Because there 
are the inertia- stored energy elements of capacitance C m and 
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FIGURE 2 | Action potential and energy function. 



inductance L m , the peak of consumption curve occurs later than 
the membrane potential. The inertia- stored energy elements in 
the biophysical model and the corresponding solving process is 
described in detail in Wang and Zhang (2007), Wang et al. (2009). 

r 0m = 0.0001 Q, ri m = 0.1 Q, r 2m = 1000 Q, r 3m = 0.1 Q, 
r m = 1000 Q, C m = 1 \iF, L m = 50 \iH, i 0m = 70.7 |jlA 

The positive and negative parts in the power curve may explain 
why cerebral blood flow is substantially increased whereas oxygen 
consumption is only slightly increased when neurons are acti- 
vated. The shape of the power curve may also explain why there 
is a synchronous effect between the external stimulus and per- 
ception generation (Igarashi et al., 2007; Wang and Zhang, 2012), 
which is still an otherwise inexplicable neurophysiological phe- 
nomenon (Dhamala et al., 2004; Haken, 2007; Wang and Zhang, 
2012). 

The decrease in energy consumption, especially the period 
of negative energy consumption, is primarily due to the local 
hyperemia induced by neural activity. Blood flow is increased 
by dilatation of blood vessels, which augments arterial blood 
flow and the local oxyhemoglobin content of brain tissue. 
Neurons absorb oxygen from oxyhemoglobin, but the oxygen 
consumption does not rise proportionately with this increased 
blood flow and oxygen supply on the millisecond timescale. 
Indeed, Fox et al. (Lin et al., 2010) observed that the oxygen 



extraction fraction (OEF) decreased from a resting value of 
approximately 40% to approximately 20% during activity asso- 
ciated with task performance. That is, the majority of oxygen 
carried in oxyhemoglobin is not metabolized by neurons, and 
only a small number of oxygen is absorbed in neurons. So the 
overall result is that neurons first exploit energy from the blood 
after the stimulus, which manifests the negative consumption 
phase. 

In terms of the relationship between energy and cognitive 
behavior, one important example is Chengyu Li's research (Li 
et al., 2009). This paper points out that cluster firing of neurons in 
cerebral cortex can shift two behavioral status of sleep and wake- 
fulness. Subthreshold and sup threshold firing states are closely 
related to neural energy. We need to stress that the energy analysis 
of structural neural networks of this paper is the basis of further 
analysis for functional neural networks, which closely related to 
higher cognitive function. 

STRUCTURAL NEURONAL NETWORKS MODEL 

The connection structure of a cortical neural network is shown 
as Figure 3, where a neuron couples with all other neurons in 
the model. The characteristics of each neuron are represented by 
the biophysical model (Figure 1). Therefore, the following neu- 
ral network structure is strictly defined on a neurobiological basis 
(Wang and Zhang, 2006, 2007). 

The connections among cortical neurons are complicated, 
one neuron couples with an estimated 10 4 peripheral neurons 
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FIGURE 3 | Neural network connection structure. 



(Huang, 2000; Singer, 2009). The global- connected neural net- 
work is composed of 15 neurons in Figure 3 in order to under- 
stand the energy- coding mode in the stimulated neural network. 
The neurons are coupled with bidirectional asymmetry cou- 
pling strengths. For example, the 1th neuron is coupled to the 
2th neuron with coupling strength 0.15, and the 2th neuron is 
coupled to 1th neuron with coupling strength 0.22. A synaptic 
coupling strength between neurons is uniformly distributed in 
the statistical sense (Rubinov et al, 2011), and we assume that 
a synaptic coupling strength is uniformly distributed in [0.1, 0.3]. 
Let the coupling strength matrix: 



Wi, l Wi, 2 • • • W\,n 
w _ W 2 ,i ' • W 2 , n 

Wn, 1 Wn,2 w n,n 



Wjj is coupling strength when the ith is coupled to jth. 
The network operates as follows: 



Sj(t) = W xQ(t- xf 

Im,j(t) = 



n [km (j - l) Sin(0) m (/ - 1) 
iml + E ~ tj-i))] (?) 

i= 1 +iom(n) sin(co OT (n)(f - t„)) if Sj(t) > th { 
iml if Sj(t) < th 



Substituting I m j(f) into Equation (5), we obtain Uo m j(t). 

Substituting Equation (5) into Equation (4), we obtain 
Pm,j(t)> where Sj(t) is the sum of the stimulations for the 
jth neuron at time t; Q(t — x) = [Qi(t — t), Q 2 (t — t), . . . , 



Qj(t — x), . . . Q n (t — x)] indicates neuron-firing states, which 
take values of 0 at resting and 1 at firing. 

POWER-CONSUMING PROPERTY DURING OSCILLATION IN 
VARIOUS PARAMETERS 

THE ENERGY-CONSUMING PROPERTY OF NEURAL NETWORK UNDER 
INSTANTANEOUS STIMULUS 

Given the biological neural network model of 15 neurons, the 
total energy consumed by the overall network and the firing 
record under the condition of instantaneous stimulation is shown 
in Figures 4A,B. 

In Figure 4, the 1st and 2nd neuron are instantaneously stimu- 
lated with —40 mV strength for 0.1 ms. These two neurons conse- 
quently fired under the stimulation, which caused the subsequent 
cascade reaction. Figure 4A shows the total power consumed by 
the overall network in 30-ms simulation time. Where the coupling 
strength is uniformly distributed in [0, 1]; the transform delay of 
synapses is uniformly distributed in [0.3 ms, 0.18 ms]. The neural 
impulse is record in Figure 4B, the black dot at coordinate (t, j) 
represents neuron j firing a neural impulse at the moment of t. 
Accordingly, the more streak-like the recording is, and the nar- 
rower the streak is, the higher the synchronism of the network 
behavior. 

According to Figure 4, the total energy- consuming curve 
shows no prominent periodicity over time. The negative energy 
existed at the beginning of the oscillation but was not absorbed 
subsequently. The simulation result of Figure 4 apparently shows 
that the oscillation of the network is not synchronous. In order 
to quantitatively estimate the synchronicity of the network activ- 
ity under instantaneous stimulation, the mean-max correlation 
coefficient (MCC) was employed. It is well-known that correla- 
tion coefficient is used to determine the degree of synchronicity. 
Given a biological neural network that consists of n neurons, a 
n x n correlation coefficient matrix can be calculated. Let this 
matrix be C, and the element at ith row and jth column be 
Cy, hence qj stands for the correlation coefficient of ith and 
jth neurons' membrane potentials. Previous research indicates 
that once the network synchronized, there will be at least two 
oscillating groups in the steady state. As the first neural group 
fires, other groups are stimulated. The simulation is accepted 
after a tiny interval at the presence of transform delay, then 
the stimulated groups consequently fire, while the first group 
came to rest at this period and is re-simulated by the cur- 
rently firing group. In this way, a coupling oscillation is formed. 
Therefore, the membrane potential of a neuron i is different 
from other n — 1 neurons, we choose the max correlation coef- 
ficient (except the correlation coefficient between the neurons 
membrane potential and itself where ca = 1) to outline the syn- 
chronization degree of the neural group containing neuron i. 
Let this the max correlation coefficient of neuron i be p^ax- 
Accordingly, an average value - X!f= l Pmax can ^ e yielded from 
P 1 o n . 

rmax' • • • ' rmax' 

The MCC is defined as follow: 



Pmean 



J2?= i max(Q 5 1, Q, i, . . . Qj, . . . , Q, „) 



N 



(i # j) 



(8) 
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FIGURE 4 1 (A) Total energy consumed by the overall network with 15 
neurons under the condition of instantaneous stimulation. (B) Overall 
network firing record with 15 neurons under the condition of 



instantaneous stimulation (the coupling strength is uniformly distributed 
in [0, 1]; the transform delay of synapsis distributed in [0.3 ms, 
0.18ms]). 



Where C h j is the correlation coefficient between the membrane 
potentials of neuron i and j, namely: 

n _ El=i\(Vi(^t)-m\Vj(kAt)-Vj\ 
jE T k =i(Vi(k&t) - Vi) 2 (Vj(kAt) - Vj) 2 

Note that the result of numerical simulation is a discrete sequence 
over time. In Equation (9), At stands for the sampling interval, 
and Vj stands for the mean value of the ith neuron's mem- 
brane potential. It is obvious that as the correlation coefficient 
between two neurons approaches 1, the stronger the degree of the 
synchronization is. 

The correlation coefficients among neurons can deal with 
both partial synchronization and coexisting synchronization in 
multiple clusters. For each neuron, we choose the maximum of 
correlation coefficients with other neurons, and then calculate 
the mean value of these maximums. The approaching to 1 of 
Pmean denotes that network's behavior is approaching the state 
of coexisting synchronization in multiple clusters; on the other 
hand, the approaching t 0 of p mea n denotes that only a subset of 
neurons are in a synchronized state, which means no coexisting 
synchronization in multiple clusters in a network's behavior. In 
Figure 4 p mea n = 0.58145, which is in accordance with the firing 
record showing that no salient synchronization could be drawn. 
Figures 5A-D shows the total energy consumed by the overall 
network in 30-ms simulation time when the neuron number is 
increased. The firing records are shown in Figures 6A-D. 

In Figure 5, the 1st and 2nd neuron are stimulated at the 
moment t = 0. The stimulation intensity is —40 mV, and the 
stimulation time is 0.1 ms, thus leading the stimulated neurons 
to fire action potentials. The energy consumed by the overall net- 
work is recorded, the coupling strength is uniformly distributed 
in [0, 1], and the transform delay of synapses is distributed in 
[0.3 ms, 0.18ms]. 



Figure 5 shows that as the number of neurons increases, the 
periodicity becomes salient, and a negative energy component 
emerges under the situation that coupling strength and transform 
delay remain steady. Meanwhile, the synchronous motion of the 
neural firing frequency tends to be prominent as the number of 
neurons increases. 

For facilitation of quantitative analysis, an index is proposed to 
outline the characteristics of energy distribution: negative energy 
ratio a(t). It is defined as the specific value in the absolute value 
of negative energy and the sum of the absolute value of negative 
and positive energy during the time period form moment 0 to 
moment t. That is, 

ct(f) = x lOOo/o (10) 

^positive 
« ft 

^negtive = / ?iW X Sgn(-Pi(t))dt (11) 

i=l Jo 

n ft 

^positive = 

J2 / Pi(t) x sgn(Pi(t))dt (12) 
i=i Jo 

Where P;(f) is the power consumed by neuron i at t moment, and 
the integration of P{(f) in [0, t] stands for the energy consumed in 
this period, sgn(-) is signal function, which is defined as sgn(x) = 

{ X > ^ . a(t) is a foundation of neural energy analysis. The 
0, x < 0 w by y 

corresponding negative ratio and MCC for Figures 5A,B is shown 
in Table 1. 

According to the simulation results of Figures 5, 6 and Table 1, 
the firing recording tends to be streak-like as the number of 
neurons increases, and greater streak width indicates that the 
synchronization tends to strengthen the oscillation. The numer- 
ical calculation in Figures 5A-D shows that the correspond- 
ing MCC efficiency increases. The negative energy ratio is also 
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FIGURE 5 | (A) The network consists of 30 neurons, pmean(a) = 0.76423. of 100 neurons, pmean(c) = 0.9476. (D) The network consists of 500 

(B) The network consists of 50 neurons, pmean(b) = 0.8555 (the neurons, pmean(d) = 0.9998 (the coupling strength is uniformly 

coupling strength is uniformly distributed in [0, 1]; the transform delay distributed in [0, 1]; the transform delay of synapsis distributed in 

of synapsis distributed in [0.3 ms, 0.18 ms]). (C) The network consists [0.3 ms, 0.18 ms]). 



monotonously increased with the number of neurons. In addi- 
tion, the energy distribution outlined by the power curve is closely 
related to network synchronization. Figures 5, 6 demonstrate 
that this correlation is monotonic. That is to say, in a coupled 
neuron network, a larger number of neurons makes it easier 
for the network to synchronously oscillate, and, simultaneously, 
the negative energy ratio increases (energy is the integration of 
power, in Figure 5, it is the area cut by power curve and the 
x-axis). 

Figures 5D, 6D show 500 neurons in synchronous oscillation; 
the periodicity of power consumed is remarkably salient in syn- 
chronous oscillation, and the negative energy power increases 
prominently. The streak of neural firing recording is slender and 
dark feature. The period form moment 0 to 14 ms is a tran- 
sient process: the peak of power was lower, and the negative 
energy augmented this feature. After 14 ms, the energy distribu- 
tion presented prominent periodicity, suggesting that the network 
achieved a synchronous oscillation state. 



CORRELATION BETWEEN COUPLING STRENGTH AND ENERGY 
DISTRIBUTION 

Given a neural network comprised of 200 neurons, the transform 
delay is uniformly distributed in [0.3 ms, 1.8 ms]. Figures 7A-D 
shows the total power consumed by the overall network when the 
coupling strength is distributed in [0, 0.01], [0, 0.1], [0, 0.5], and 
[0, 1], respectively. The firing record during oscillation is shown 
in Figures 8A-D. 

In the simulation of Figures 7, 8, the 1st and 2nd neurons 
are stimulated at the moment t = 0, the stimulating intensity 
is — 40 mV, and the stimulation lasted 0.1 ms. In order to avoid 
overall tranquillization (where the coupling strength between the 
stimulated neurons and the others is so slight that other neu- 
rons cannot be stimulated to fired), the coupling strength between 
the stimulated neurons (1st and 2nd neuron) maintained a uni- 
form distribution in [0, 1]. At the same time, coupling strength 
among other neurons is uniformly distributed in [0, 0.01], [0, 
0.1], [0, 0.5], and [0, 1], respectively, for Figures 7A-D. The 
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FIGURE 6 1 (A) The firing record for network with 30 neurons. 
(B) The firing record for network with 50 neurons (the coupling 
strength is uniformly distributed in [0, 1]; the transform delay of 
synapsis distributed in [0.3 ms, 0.18 ms]). (C) The firing record for 



network with 100 neurons. (D) The firing record for network 

with 500 neurons (the coupling strength is uniformly distributed 

in [0, 1]; the transform delay of synapsis distributed in [0.3 ms, 
0.18ms]). 



Table 1 | The corresponding negative ratio and mean-max correlation 
coefficient for Figures 5A,B. 



Neuron amount 


30 


50 


100 


500 


a(f) (%) 


0.1526 


0.2362 


0.4532 


1 .9457 


Pmean 


0.7642 


0.8556 


0.9746 


0.9998 



The distribution interval of couple strength: [0, 1] 

The distribution interval of transform delay [0.3ms, 1.8ms]. 



corresponding negative energy ratio and MCC are displayed in 
Table 2. 

According to the calculations of Figures 7, 8, and Table 2, as 
the coupling strength increases, the periodicity, negative energy 
ratio, and MCC of the total power curve and the synchronization 
reflected by firing recordings exhibited a consonant monotonic- 
ity. In other words, the total power curve usually appears to be 
approximately periodic when the oscillations reach a steady state, 



and the larger values of negative energy ratio and MCC and the 
salience of synchronization are reflected by the firing recordings. 

CORRELATION BETWEEN EXCITABLE SIGNAL TRANSFORM DELAY AND 
ENERGY DISTRIBUTION FEATURE 

Given the number of neurons and the constant distribution 
interval of coupling strength, the feature of total energy con- 
sumed by the overall network varies as the signal transform delay 
changes. Figures 9A-D shows that the total power consumed by 
the overall network as the distribution interval of signal varies 
under the condition that the number of neurons is maintained 
at 300, and the coupling strength is uniformly distributed in 
[0, 1]. Figures 10 A-D shows the corresponding firing impulse 
recordings during the oscillation. 

In the simulation of Figures 9, 10, the 1st and 2nd neuron are 
stimulated at the moment t = 0 at — 40 mV, and the stimulation 
lasted 0.1 ms. The coupling strength is uniformly distributed in 
[0, 1] while the signal transform delay is uniformly distributed in 
[0.1ms, 1.6ms], [0.2ms, 1.7ms], [0.3ms, 1.8ms], and [0.5ms, 
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FIGURE 7 | (A) Total power consumed by the overall network the coupling 
strength is uniformly distributed in [0, 0.01] pmean = 0.8567. (B) Total power 
consumed by the overall network the coupling strength is uniformly 
distributed in [0, 0.1] pmean = 0.8567. The amount of neuron is 200, and the 
transform delay is uniformly distributed in [0.3 ms, 1.8 ms]. (C) Total power 



10 15 

time/ms 



consumed by the overall network the coupling strength is uniformly 
distributed in [0, 0.05] pmean = 0.9470. (D) Total power consumed by the 
overall network the coupling strength is uniformly distributed in [0, 1] 
p mean = 0.9827. The amount of neuron is 200, and the transform delay is 
uniformly distributed in [0.3 ms, 1.8 ms]. 



2 ms], respectively. The corresponding negative energy ratio and 
MCC are displayed in Table 3. 

According to Figure 9 and Table 3, the negative energy 
ratio increases as the distribution interval of transform delay 
approaches zero. The direction of the correlation coefficient, how- 
ever, cannot guide any persuasive conclusions. This is because the 
MCC merely reflects the local synchronization among one spe- 
cific neural population. As multi-population oscillation is formed 
under the instantaneous stimulation, the intergroup synchroniza- 
tion cannot be reflected by MCC. According to the neural impulse 
recorded in Figure 10, it is obvious to note that the intergroup 
synchronization decreases as the distribution interval of trans- 
form delay is further from zero. In Figure 10C, it can be easily 
perceived that in two stripes with a relative wider interval, called a 
first- order stripe, exist two stripes with a relative narrow interval, 
called a second- order stripe, which means that two "micro" oscil- 
lation groups exist in each "macro" oscillation group. Because the 



synchronous degree between the micro oscillation group corre- 
spondence to the second-order stripe is high, no definitive result 
can be achieved by quantifying the correlation coefficient. With 
regard to thee measure of neural energy, however, the intergroup 
non- synchronization can be effectively clarified by the negative 
energy ratio. 

The reason why energy analysis can quantitatively analyze 
the intergroup synchronization is that the negative power only 
emerges during a tiny period at the beginning of the neural 
action process, and it takes a trivial ratio of the total energy. 
If the power of neurons is linearly superimposed, the nega- 
tive component will be offset by the positive component when 
the oscillation synchronization is not salient, which leads to no 
negative component in the total power curve. Moreover, as the 
intergroup oscillation is formed in the network under instan- 
taneous stimulation, offsetting occurs if the phase difference is 
distinct. 
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FIGURE 8 | (A) The firing impulse record of the network's neural activity impulse record of the network's neural activity the coupling strength is 
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impulse record of the network's neural activity. The coupling strength is network's neural activity the coupling strength is uniformly distributed in 
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THE RELATIONSHIP BETWEEN NEURON NUMBER AND 
NETWORK ENERGY FEATURES 

THE RELATIONSHIP BETWEEN THE NEURON NUMBER AND THE 
ENERGY FEATURE 

For the study of the relationship between the neuron number and 
the energy feature, let the signal transform delay be uniformly 
distributed in [0.2 ms, 1.7 ms] (Dhamala et al., 2004) and the cou- 
pling strength is uniformly distributed in [0, 0.8]. The negative 
energy ratio a(t) and MCC are calculated as the neuron number 
varies from 50 to 1000, these are shown in Figures 11A,B- 
Each point in Figure 1 1 was obtained as follows: 

(1) Let signal transform delay uniformly distribute in [0.2 ms, 
1.7 ms] and the coupling strength uniformly distribute in [0, 
0.8]; 

(2) Stimulate any two neurons in the network at —40 mV for 
0.1 ms; 

(3) Calculate the total power consumed by the overall neuron 
population P m (t); 



(4) Obtain a and p mea n by the way of Equations (8)— (12); 

(5) Repeat the process above 10 times and yield 10 of a and p mea n 
to obtain the average of these 10 a and pmean? termed a and 

Pmean 

(6) a and p mean are the values of vertical coordinates in 
Figures 10A,B. 

According to Figure 11, the relationship between the neuron 
number and the negative energy is monotonic, which means 
as the number of neurons increases, the negative energy ratio 
also grows. When there are more than 400 neurons, "saturation" 
occurs, which is a consequence of the saturation of synchro- 
nization. In the mean time, the MCC showed a similar feature 
with the variation of the neuron amount. This indicates that the 
energy-based neural coding theory is also capable of represent- 
ing network activity synchronization and is highly accordant with 
the traditional measure of correlation coefficient. In fact, energy- 
based neural coding is superior to correlation coefficient-base 
analysis. On one hand, the MCC curve comes to saturation 
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Table 2 | The corresponding negative energy and mean-max 
correlation coefficient. 



Coupling strength 


[0, 0.01] 


[0,0.1] [0,0.5] 


[0,1] 


a(f) (%) 
Pmean 


6.5122 x 1CT 3 
0.8567 


1.0094 x 10" 2 0.4919 
0.9332 0.9470 


1 .4588 
0.9827 


The amount of neurons: 200. 
Distribution Interval (uniform) [0.3 ms, 1. 


8 ms]. 




Table 3 | The corresponding negative energy and mean-max 
correlation coefficient. 


The distribution 
interval of 
transform delay 


[0.1 ms, 
1.6 ms] 


[0.2 ms, [0.3 ms, 
1.7 ms] 1.8 ms] 


[0.5 ms, 
2 ms] 


a(f) (%) 
Pmean 


2.731 
0.9980 


1 .7609 1 .6832 
0.9982 0.9986 


0.6131 
0.9671 



The amount of neurons: 300. 



The distribution interval of 10, 1]. 

when the neurons outnumber 200, whereas saturation occurred 
as the number of neurons exceeded 400 for the negative energy 
ratio. Moreover, the negative energy ratio maintained a slope in 
the saturation area. On the other hand, the superiority can be 
employed to study the influence of multiple variable parameters 
on synchronization, as Figure 12 shows. 

In Figure 12, the signal transform delay is uniformly dis- 
tributed in [0.1ms, 1.6 ms] (black dot), [0.3 ms, 1.8 ms] (blue 
dot), and [0.5 ms, 2 ms] (green dot), and each point (both a(t) 
and pmean) is achieved in the same way as for Figure 11. 

According to Figure 12, the negative energy ratio and MCC 
both reflected the differences of the neuron number and sig- 
nal transform delay. When there were more than 200 neurons, 
however, the MCC curves with transform delay intervals of 
[0.1ms, 1.6 ms] and [0.3 ms, 1.8 ms] are almost indistinguish- 
able. Moreover, the MCC curve with transform delay uniformly 
distributing in [0.5 ms, 2 ms] is relatively close to the other two 
curves, which is a poor distinguishing feature for the variation 
of signal transform delay. In this sense, the correlation coefficient 
can merely reflect the variation for single parameter. Nevertheless, 
negative energy ratio can distinctly and simultaneously reflect the 
signal transform delay and the neuron number; the three curves 
did not converge. It can be easily noted that the negative energy 
power decreased as the distribution interval of signal transform 
delay moved away from zero. 

According to the discussion above, the measure of energy 
can yield more valid, valuable conclusions dealing with mutually 
coupled neuron populations' synchronization and the parameter 
distribution. Moreover, it is capable of addressing many complex 
problems that traditional methods cannot. 

THE RELATIONSHIP BETWEEN NEURAL ENERGY AND THE COUPLING 
STRENGTH 

Given a neural network with 200 neurons and the signal trans- 
form delay uniformly distributed in [0.3 ms, 1.8 ms], we studied 



the relationship between neural energy feature and the cou- 
pling strength. The total power consumed by the overall network 
was numerically calculated with varying coupling strength dis- 
tributing intervals, and the negative ratio a(t) was obtained as 
well as the MCCs pmean- The curve of a(t) and p mea n over the 
variation of coupling strength distributing interval is plotted in 
Figures 13A,B. 

Each point in Figure 13 is obtained as follows: 

(1) Let the coupling strength among neurons uniformly dis- 
tribute in [0, 0.1 + A], where A is 0, 0.01, 0.02,..., 1, 
totally 100 value, and each A correspond to a point in the 
figure. There is respectively 100 points in Figures 13A,B. 

(2) Employ a neural network with 200 neurons and a signal 
transform delay uniformly distributed in [0.3 ms, 1.8 ms]; 

(3) Stimulate any of two of the neurons in the network 
at -40 mV for 0.1ms; 

(4) Calculate the total power consumed by the neuron popula- 
tion P m (t); 

(5) Calculate a and pmean in the way of Equations (8)— (12) 

1. Repeat the process above 10 times to yield 10 a and pmean and 
obtain the averages: a and p mea n- 

2. a and p mean are the vertical coordinate values in 
Figures 13A,B. 

According to Figure 13, both the negative energy ratio and the 
MCC are monotonic with the variation of the coupling strength: 
the bigger value of the coupling strength is distributed, the big- 
ger values of the negative ratio and the MCC yield. However, the 
changing of slopes is different: it is clear that the curve of nega- 
tive energy ratio has two inflection points. Based on the analysis 
above, it is result in the asynchronous intergroup oscillation. That 
means, the energy-based neural coding reflects not only the syn- 
chronizing affection caused by coupling strength, but the change 
of the synchronization of intergroup oscillation as the coupling 
strength varies. 

THE RELATIONSHIP BETWEEN THE SIGNAL TRANSFORM DELAY AND 
THE ENERGY FEATURE 

In order to study the relationship between energy feature and 
the signal transform delay, a network with 300 neurons was 
employed, where the coupling strength was distributed in [0, 
1]. The total power consumed was calculated under the instan- 
taneous stimulation with various distributing interval of signal 
transform delay, and therefore the curve of negative energy ratio 
a(f), shown in Figure 14A was obtained. 

Each point in Figure 14 was obtained as follows: 

(1) The signal transform delay is distributed in the interval set 
[0.01 + A, 1.5 + A], where A is 0, 0.01, 0.02,..., 0.7, 
respectively, for each interval and corresponds to one dot in 
Figure 14. There are a total of 70 dots in Figures 14A,B. 

(2) Let the neuron number be 300 and the coupling strength 
uniformly distributed in [0, 1]. 

(3) Stimulate any two neurons in the network at — 40 mV for 
0.1ms. 
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FIGURE 9 | (A) Signal transform delay is uniformly distributed in [0.1 ms, 
1.6 ms] pmean = 0.9980. (B) Signal transform delay is uniformly distributed 
in [0.2 ms, 1.7 ms] pmean = 0.9982. The amount of neuron is 300, and the 
coupling strength is uniformly distribute in [0, 1]. (C) Signal transform delay 



10 15 

time/ms 



is uniformly distributed in [0.3 ms, 1.8 ms] pmean = 0.9980. (D) Signal 
transform delay is uniformly distributed in [0.5ms, 2 ms] pmean = 0.9982. 
The amount of neuron is 300, and the coupling strength is uniformly 
distribute in [0, 1]. 



(4) Calculate the total power consumed by the neuron popula- 
tion P m (t); 

(5) Calculate a and p mea n in the way of Equations (8)— (12) 

(6) Repeat the process above 10 times to yield 10 a and pmean and 
obtain the averages: a and p mean - 

(7) a and p mean are the vertical coordinate values in 
Figures 13A,B. 



In Figure 14, the curve of the negative energy ratio and the MCC 
can be divided into four sections. The first corresponds to the 
horizontal axis value from 0 to 0.2. In this section, the nega- 
tive energy ratio decreased as the interval of signal transform 
delay moved away from zero. The MCC, however, shows no such 
trend saliently because MMC cannot reflect intergroup oscillation 
synchronization, whereas the negative energy ratio can. The sec- 
ond section corresponds to the curve where the horizontal axis 
values range from 0.2 to 0.4. In this section, neither the nega- 
tive energy ratio curve or the correlation coefficient prominently 
change as the distributing intervariation of the signal transform 
delay. It can be inferred that the signal transform delay only has 



a tiny influence on network activity synchronization. The third 
section corresponds to the horizontal axis values from 0.4 to 
0.5. In this section, the negative energy curve and the correla- 
tion coefficient both dropped drastically. In the fourth section, 
the horizontal axis value is greater than 0.5, the negative energy 
ratio maintained near zero, while the correlation coefficient kept 
a relatively large value. This phenomenon illustrates that inter- 
group synchronization is faint, but the internal synchronization 
was well- sustained. 

According to the analysis above, it can be seen that the relation- 
ship between the negative energy ratio and the signal transform 
delay is sectional. This analysis, combined with the correlation 
coefficient, shows that the synchronization of the network activity 
is discontinuously influenced by variation in the signal transform 
delay. Moreover, compared with the correlation coefficient, a tra- 
ditional measure, more detailed information can be generated 
with energetical measurement. 

ENERGY FEATURES UNDER MULTIPLE VARIABLE PARAMETERS 

The analysis of the relationship between the negative energy ratio 
and the correlation coefficient mainly concerned single variable 
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FIGURE 10 | (A) The firing impulse record of the network's neural activity 
signal transform delay is uniformly distributed in [0.1 ms, 1.6 ms|. (B) The 
firing impulse record of the network's neural activity signal transform delay is 
uniformly distributed in [0.2 ms, 1.7 ms]. The amount of neuron is 300, and 
the coupling strength is uniformly distribute in [0, 1]. (C) The firing impulse 



record of the network's neural activity signal transform delay is uniformly 
distributed in [0.3 ms, 1.8 ms]. (D) The firing impulse record of the network's 
neural activity signal transform delay is uniformly distributed in [0.5 ms, 2 ms]. 
The amount of neuron is 300, and the coupling strength is uniformly 
distribute in [0, 1]. 



parameters. The two-variable parameters circumstance is shown 
in Figures 15, 16. 

In Figure 15, the neuron number changed gradually from 50 
to 500, and the distributing interval of the coupling strength 
changed gradually from [0, 0.02] to [0, 1]. The signal trans- 
form delay was uniformly distributed in [0.3 ms, 1.8 ms]. For 
each parameter set, any two neurons are stimulated for 0.1ms 
at — 40 mV to force them fire, and we could determine the 
total power consumed in 30 ms and the membrane potential. 
Therefore, the negative energy ratio a and the correlation coef- 
ficient pmean are obtained. Each point in Figure 15 is the mean 
value of 10 simulations. 

According to the simulation results, it can be inferred that: 

(1) Increasing the neuron number and the coupling strength 
both lead to a negative energy ratio and an increasing cor- 
relation coefficient with a coherent trend. Considering the 
analysis of single -variable parameters, the effect of the neu- 
ron number and the coupling strength can be determined. 



(2) In Figure 14A the negative energy ratio gradually increased. 
However, in Figure 15 A, the correlation coefficient enhanced 
drastically at first and then came to saturation as the hor- 
izontal and ordinate value increased. In other words, com- 
pared with traditional correlation coefficients, the quantita- 
tive measure of energy method possesses better hierarchy and 
better reflects the relationship between itself and the network 
parameters. 

The conclusion above is consistent with the previous discussion. 
Multiple oscillation groups will be emerge under the instanta- 
neous stimulation after the network achieves a stable state, and 
the intergroup oscillation is not always synchronous and is sensi- 
tive to the distributing interval parameter. It is hard to determine 
if this problem can be analyzed by the traditional correlation coef- 
ficient measure. However, with the energy method studied in this 
paper, a better effect can be achieved. 

In Figures 15, 16, the distributing interval of the coupling 
strength changed gradually from [0, 0.02] to [0, 1], and the 
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distributing interval of the coupling strength changed grad- is simulated, as well as the membrane potential. Therefore, 

ually from [0.01, 0.01 +1.5] to [0.5, 0.5 +1.5]. The num- the negative energy ratio a and the correlation coefficient 

ber of neurons was maintained at 200. For each parameter pmean are obtained. Each point in Figure 15 is a mean value 

set, any two neurons are stimulated for 0.1ms at — 40 mV to of 10 simulations. It can be seen that a similar result was 

force them to fire; hence, the total power consumed in 30 ms achieved. 
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FIGURE 13 | (A) The curve of a(f) over the variation of coupling strength distributing interval. (B) The curve of pmean over the variation of coupling strength 
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In Figure 17, the distributing interval of the coupling strength 
changed gradually from [0, 0.02] to [0, 1] and the neuron num- 
ber changed gradually from 50 to 500. The coupling strength was 
uniformly distributing in [0, 1]. For each parameter set, any two 
neurons are stimulated for 0.1 ms at — 40 mV to force them fire, 
and the total power consumed in 30 ms and the membrane poten- 
tial are simulated. Therefore, the negative energy ratio a and the 
correlation coefficient p mea n are obtained. Each point in Figure 15 
is the mean value of 10 simulations. 

With the calculation in Figure 17, the discussion about sig- 
nal transform delay can be further explained. Figure 17B shows 
that signal transform delay has a small influence on MCC but 
has a distinct effect on the negative energy ratio. This is because 
the negative energy ratio reflects the synchronization activity to 
a great extent, which the MCC is rarely capable of. Hence, we 
found that the intergroup asynchronization was mainly due to 
long signal transform delay. 

CONCLUSIONS 

In this paper, we studied the relationships among neural energy 
features, network parameters, and oscillating characteristics. 
What should be stressed is that the value of a parameter can 
vary from one specific function to another, but the relationships 
remain coherent. Specifically, they are: 

( 1 ) The relation between negative energy ratio and the amount of 
neurons is monotonic increasing, as well as the synchronous 
oscillation (measured by correlation coefficient). 



(2) The relation between negative energy ratio and the coupling 
strength is monotonic increasing, as well as the synchronous 
oscillation (measured by correlation coefficient). 

(3) The relation between negative energy ratio and the time delay 
of transmitting signal is monotonic decreasing, as well as 
the synchronous oscillation (measured by correlation coef- 
ficient). 

(4) Negative ratio can better represents the variation of syn- 
chronous activity due to the multiply changed parameters. 

(5) Negative ratio can better represents the intergroup asynchro- 
nization. 

The study method proposed in this paper represent the net- 
works' parameters and its activity with respect to the energy 
method. It is a novel method rather than a precise conclusion. But 
still, it can quantitatively represent the status of oscillation with 
respect to synchronization or asynchronization. Multiple variable 
parameters will be an important focus of further research. 

Though the networks studied in this paper is structural, it 
is an important foundation for further development for energy 
method to address on functional networks, which is closely 
related to cognitive and behavior problems. And this is another 
focus in the future. 

The energy method is a new theory to investigate the global 
behavior of brain activity. One significant characteristic is the 
superposition property. In experimental neuroscience, record- 
ing all neural impulses, even in a small encephalic region, is 
virtually impossible. However, it is possible to measure the 
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FIGURE 15 | (A) The negative energy ratio as the neuron amount and the coupling strength simultaneously changes. (B) Correlation coefficient as the neuron 
amount and the coupling strength simultaneously changes. 



total energy of the region. Based on this, the study method 
proposed in this paper is a effective tool for experimen- 
tal neuroscience; parameter distribution of a given encephalic 
region can be estimated by the energy method proposed 



in this paper. This method can also be used to study dif- 
ferences between encephalic regions. It is of more practical 
value and promotive value compared with traditional numerical 
analysis. 
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FIGURE 16 | (A) The negative energy ratio as the coupling strength and the signal transform delay simultaneously changes. (B) The correlation coefficient as 
the coupling strength and the signal transform delay simultaneously changes. 
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FIGURE 17 | (A) The negative energy ratio as the coupling strength and the signal transform delay simultaneously changes. (B) The correlation coefficient as 
the coupling strength and the signal transform delay simultaneously changes. 



The theoretical study of nervous system coding and decoding 
is in bottleneck period, and new methods are needed to augment 
traditional ones and achieve more valuable results. It is important 
to mention that the energy coding theory is in its infancy, and 
there are numerous problems that remain to be solved. However, 
the energy method has significant potential and could be useful 
for addressing questions that traditional neuroscientific methods 
have been unable to elucidate. 
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